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Abstract 

Estimates for electroproduction of pseudoscalar mesons at small 
values of skewness are presented. Cross sections and asymmetries for 
these processes are calculated within the handbag approach which is 
based on factorization in hard parton subprocesses and soft generalized 
parton distributions (GPDs). The latter are constructed from double 
distributions. Transversity GPDs are taken into account; they are 
accompanied by twist-3 meson wave functions. For most pseudoscalar- 
meson channels a combination of Ht and Et plays a particularly 
prominent role. This combination of GPDs which we constrain by 
moments obtained from lattice QCD, leads with the exception of the 
7r + and n' channels, to large transverse cross sections. 
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1 Introduction 

In a recent paper [1] we have investigated hard exclusive electroproduction 
of positively charged pions within the framework of the handbag approach in 
the kinematical range of small skewness, £, and small invariant momentum 
transfer, —t but large photon virtualities, Q 2 . As in our previous studies of 
vector-meson electroproduction [2, 3] the partonic subprocess is calculated 
within the modified perturbative approach [4] in which quark transverse de- 
grees of freedom as well as Sudakov suppressions are taken into account. 
In other words the transverse size of the produced meson is not ignored as 
in the collinear (leading-twist) approach. On the other hand, the partons 
entering the subprocess are viewed as being emitted and reabsorbed by the 
nucleon collinearly to the nucleon momenta. The GPDs which embody the 
soft physics, are constructed with the help of double distributions. In [1] the 
sketched approach has been used to analyze the HERMES data on n + elec- 
troproduction [5, 6, 7]. Besides the use of the modified perturbative approach 
this analysis differs in the following aspects from previous studies based on 
the collinear approximation (e.g. [8] - [11]): 

• The full electromagnetic form factor of the pion, F 7T {Q 2 ), is taken into 
account (see also Ref. [12]) as it is measured by the F w collaboration 
[13] at Jefferson Lab in just the same process 

• Besides the pion pole there is an extra contribution to the GPD E. 

• There are substantial contributions from transversely polarized pho- 
tons. This is particularly obvious from the sin s -moment, A s ^ a , of 
the 7r + cross section measured with a transversely polarized target [6]. 
As it is argued in [1] the contributions from transversely polarized pho- 
tons can be calculated within the handbag approach as a twist-3 ef- 
fect consisting of a twist-3 pion wave function and the leading-twist 
transversity GPD H T . 

Here, in this work we are going to extend the analysis performed in [1] to 
other pseudoscalar meson channels, namely 7r°, K and 77(77') production. For 
these processes there are no small-skewness data available as yet but they 
may be measured at the upgraded Jefferson Lab facility or by the Compass 
experiment. In so far we believe that estimates for various pseudoscalar 
meson channels are of interest and timely. While for 7T + production only 
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the isovector combinations of the GPDs, = F u - F d (F = H,E,H T ), 
contribute, different flavor combinations are relevant for the other channels. 
This necessitates a careful reexamination of the parameterizations of the 
GPDs with particular regard to those of the individual flavors. In some 
cases, as for instance for E, we have to revise slightly the parameterizations 
proposed in [1]. In a recent lattice-QCD study [14] large moments of the 
transversity GPD combination 

Ex — 2Hj< + Ex (1) 

has been found. With regard to this result we will also examine the role 
of this GPD which has not been taken into account in [1], in order to find 
out whether or not it provides substantial effects in observables for meson 
electroproduction. 

The plan of the paper is the following: In the next section we will sketch 
the handbag approach to meson electroproduction including twist-3 effects. 
In Sect. 3 we are going to present the parameterizations of the GPDs and 
compare them with recent results from lattice QCD in detail. Predictions for 
electroproduction of pseudoscalar mesons will be presented in the subsequent 
sections, for pions in Sec. 4, for rj and rj' in Sect. 5 and for kaons in Sect. 6. 
In Sect. 7 the asymmetries obtained with either a longitudinally polarized 
beam or target will be discussed. The paper is closed with a summary (Sect. 
8). 



2 An outline of the handbag approach 



For details of the approach we are going to use it is referred to our previous 
papers, e.g. [1, 2, 3]. Here, we only sketch the basic facts. It is to be stressed 
that we consider the kinematical region of small £ and small —t but large 
Q 2 and large photon-proton c.m.s. energy, W. Terms of order {\f—i/Q) n 
(n > 1) are neglected throughout. 

The helicity amplitudes for electroproduction of pseudoscalar mesons, 
7*p — > PB, through longitudinally polarized photons read 
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Helicities are labeled by their signs or by zero. The usual abbreviation t' = 
t — t is employed where [11] 

_ (m 2 + M 2 )£ 2 + (M 2 - m 2 )£ 
*o - ~ z j^7|2 w 

is the minimal value of — t corresponding to forward scattering. The mass of 
the nucleon (meson, final state baryon) is denoted by m(mp, M), eo is the 
positron charge and the skewness is related to Bjorken-x by 

XB -Ml + m 2 P /Q 2 ] . (4) 



The pole contribution, see Fig. 1, occuring in n + and K + production has 
the residue 

Pp = gp P BF PpB (t)Fp(Q 2 ) , (5) 

where gppp,Fp P B{t) is the coupling of the meson to the proton-baryon vertex 
and Fp represents the electromagnetic form factor of the meson for which 
we use the experimental values. For tt° and r\ production pp is zero. Last 
not least the item {F p ) in (2) denotes a convolution of the GPD F with an 
appropriate subprocess amplitude to be calculated from a set of Feynman 
graphs of which a typical leading-order example is shown in Fig. 1, 

(^ P ) = E/ 1 dxH^ 0X (x,Z,Q 2 ,t = 0)F p (x,Z,t). (6) 

The label A refers to the unobserved helicities of the partons participating in 
the subprocess. Within the modified perturbative approach the subprocess 
amplitude reads 

Ho X ,ox = J dTd%^p(T,-h,p F )f p XfiX (x,£,T,Q 2 ,h,fi R ) 

xa s (fi R )exp[-S(T,b,Q 2 ,fi F ,fi R )} (7) 

in the impact parameter (b) space. For the Sudakov factor S, the choice of 
the renormalization {hr) and factorization (/ip) scales as well as the hard 
scattering kernels J 7 or their respective Fourier transforms F, we refer to 
Ref. [3]. The last item in (7) to be explained is \&p(r, — b) which represents 
the Fourier transform of the momentum-space light-cone wave function for 
the meson (r is the momentum fraction of the quark that enters the meson, 
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Figure 1: The pole contribution to the process 7*p — > P + B(\eft) and the 
handbag graph (right) for electroproduction of pseudoscalar mesons. Helici- 
ties are specified for the amplitude .Mo-,++- 



defined with respect to the meson momentum, the antiquark carries the 
fraction 1 — r). 

As forced by the small —t' behavior of the HERMES data [6] on A S j^ s for 
7r + electroproduction one also needs the amplitude A^o-.++- As is obvious 
from the parton helicities specified in the Feynman graph depicted in Fig. 1, 
the usual GPDs H and E which parameterize the nucleon matrix element 
for a situation where the emitted and reabsorbed partons have the same 
helicity, do not provide a contribution to .Mo-,++ with the required behavior 
.Mo-,++ — > const, for if — > 0. Angular momentum conservation forces this 
contribution to vanish in the forward direction. What is required, as has been 
discussed in detail in [1] , is a contribution from the helicity-flip or transversity 
GPD Ht [15, 16] in combination with a twist-3 meson wave function 



M 



e,Jl~e £ dx H^ 3 H£ . (8) 



This amplitude is parametrically suppressed by fip/Q with respect to the 
asymptotically leading amplitudes for longitudinally polarized photons. The 
parameter fip is large since it is enhanced by the chiral condensate 

HP= r~ , (9) 

by means of the divergency of the axial- vector current. The m 8 are the 
current-quark masses of the meson's valence quarks. The parameter \xp is 
scale dependent and evolves with the anomalous dimensions 4//5q = 12/25 
for four flavors. For the rj and rf the situation is a bit more complicated. 
Decomposing these mesons into flavor-octet and single states, one has for the 
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octet case 

2 

TYl 

H = 3 m - . (10) 

m u + m d + Am s 

The flavor-singlet 77 is not related to the chiral condensate. In this case i^rji 
is just the 771 mass which we approximate by that of the rj'. For pions, kaons 
and rjs we take a value of 2 GeV for the parameter /j, p at a scale of 2 GeV. 

The twist-3 subprocess amplitude is given explicitly in [1] for 7r + electro- 
production. Its generalization to other pseudoscalar mesons is straightfor- 
ward. 

As we mentioned in the introduction we will also examine the role of the 
GPD (1) which contributes to the amplitude 

where \x = ±1 indicates a transverse photon helicity. For our kinematical 
range of small —t' and small skewness contributions from the other transver- 
sity GPDs to the amplitudes (8) and (11) can be neglected. Note that E T is 
an odd function of £ as a consequence of time-reversal invariance [16]. The 
double-flip amplitude Aio--+ is also suppressed in the kinematical range of 
interest and neglected. 

The pion-pole contributions to the amplitudes for transversely polarized 
photons are also taken into account. Explicit expressions for these contribu- 
tions can be found in [1]. 



3 The parameterization of the GPDs 

According to the discussion presented in Sect. 2 we need to model the GPDs 
H, E, H T and E T . It is to be emphasized that our parameterizations of 
the various GPDs are optimized for small skewness and small —t. We as- 
sume flavor-symmetric sea GPDs throughout. Hence, we have to model only 
the valence-quarks GPDs since the sea contributions cancel for 7r + and K + 
production on this assumption. For the cases of tt° and 77 the sea does any- 
way not contribute. At present there is no experimental information from 
hard exclusive processes available which would allow for an examination of 
flavor-symmetry breaking effects in the sea. In recent studies of the polarized 
parton distributions [17, 18] however flavor-symmetry breaking effects in the 
sea have been found. Thus, at least for H an estimate of the sea contribution 
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on the basis of the double distribution ansatz, is possible. We found its sea 
contribution to the n + production cross section to amount to about 1%. It 
therefore can safely be neglected. 

The valence- quarks GPDs are constructed with the help of their double 
distribution representation. For the latter the familiar ansatz [19] 

f?(p, rj, t) = exp \{h - a> \np)t] F?(p, £ = t = 0) \ [(1 ~ P ^~/ ] 9(p) 

(12) 

is made which consists of the forward limit, F i: of the relevant GPD, a weight 
function and an exponential in t with a profile function parameterized in a 
Regge-like fashion, i.e. with a slope of a Regge trajectory, a[, and a slope, bi, 
of the residue function. The t — part of the Regge term is absorbed in the 
forward limit of the GPD which is given either by the polarized (Aq a (x)) or 
the transversity (5 a (:r))j3arton distribution functions (PDFs) for quarks of 
flavor a for the GPDs H a or H£. For the remaining two GPDs, E a and E?, 
the forward limit is parameterized in a fashion analogously to the PDFs 

F?(p, £ = t = 0) = N9p-^°) (1 - p)tn (13) 

with the parameters to be adjusted appropriately. 

The full GPDs are obtained from the double distributions by the integral 
[20, 21] 

f i r i-\ P \ 

F t a (x, e, t) = / dp d V 5(p + £v-x) ft (p, V , t) . (14) 

J-i J-i+\p\ 

Evolution is neglected throughout since we are only interested in scales in 
the proximity of 2 GeV which is the scale at which we quote the parameter- 
izations of the GPDs. Below we will occasionally refer to moments of GPDs, 
in particular to moments at zero skewness. They are defined by 

F in0 (t) = £ i dxx n - 1 F i (x,S = 0,t). (15) 

Before we turn to the discussion of the individual GPDs we want to 
comment on the Regge trajectories, aii(t), used by us. They are to be under- 
stood as effective trajectories describing unspecified superpositions of poles 
and cuts. They are not related to the hadronic spectrum as for instance the 
prominent p and oj Regge trajectories which contribute to the GPDs H and 
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E. The appearance of cuts in soft reactions dominated by unnatural-parity 
exchanges is a well-known fact. For definiteness let us consider the closely re- 
lated reaction photoproduction of mesons (M). For single-particle exchanges 
the corresponding helicity amplitudes satisfy the symmetry relation 3 

M V,-,. = ^m(-I)^' M^ u , (16) 

where, for vector (pseudoscalar) mesons, i] M = +(— )1 and, for a (un)natural 
parity exchange, k = +(— )1. Eq. (16) relates for instance the amplitudes 
A4q-,-+ and .Mo-,++- The first one is a double-flip amplitude which, by 
angular momentum conservation, is forced to vanish oc t' for forward scat- 
tering. Therefore, A^o-,++ also vanishes oc t' for dynamical reasons although 
it is a non-flip amplitude. Hence, for n + photoproduction, all single-particle 
exchanges vanish for forward scattering and would produce a forward dip in 
the differential cross section. This is in sharp contrast to experiment [22] - 
the 7r + photoproduction cross section exhibits a pronounced forward spike. 

3.1 H 

The forward limit of H, the polarized parton distribution 4 is taken from [23]. 
The Regge parameters for H are quoted in Tab. 1. In [24] a profile function 
for H has been proposed that is more complicated than bfj — ag, In p in (12). 
It reads 

-ofe(l - pf hip + b a 6 (l - pf + A\p(\ - pf (17) 

In [24] the parameters appearing in (17), have been determined from the 
data on the axial-vector form factor [25] through the sum rule for H: 

a' 6 = 0.9 GeV" 2 , b\ = b% = 0.59 GcV~ 2 

A\ = 1.22 GeV~ 2 A% = 2.59 GeV" 2 . (18) 

This profile function is also valid at large —t. As has been noted in [24] there 
is a strong correlation between p and t in the zero-skewness GPDs: at large 

3 In [1, 3] it is shown that this relation also holds for the handbag amplitude at the twist- 
2 level (with or without power corrections like effects and to any order of perturbative 
QCD). The contributions from the GPDs H and E {H and E) corresponds to (un)natural 
parity exchange. Contributions from Et behaves as natural parity exchanges while those 
from Ht do not possess the property (16). 

4 A suitable parameterization of these PDFs can be found in [3]. As compared to that 
work the parameterization of Au v is changed somewhat. 
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GPD 


a(0) 


a'[ GeV- 2 ] 


6[ GeV" 2 ] 


N u 


N d 


H 


0.48 


0.90 


0.59 






E n .p. 


0.48 


0.45 


0.9 


14.0 


4.0 




-0.02 


0.45 





0.78 


-1.01 


Et 


0.3 


0.45 


0.5 


6.83 


5.05 



Table 1: Regge parameters and normalizations of the GPDs, quoted at a 
scale of 2 GeV. 

(small) —t the large (small) p-region dominates. With regard to this fact 
one can view the profile function in (12) as the small —t approximation of 
the more complicated profile function (17). Indeed both the profile functions 
for H lead to practically the same results for the axial-vector form factor 
and for pion electroproduction at small —t. In Fig. 2 the first moment of 
the isovector combination of H which represents the axial form factor of the 
nucleon 

F A (t) = J^dxHQ\x,Z,t), (19) 

is shown and compared to the experimental data [25] and to results obtained 
from a recent lattice-QCD simulation [26]. Agreement between the moments 
evaluated from H and experiment is to be observed while the lattice result 
exhibits a flatter t dependence although it corresponds with the others at 
small —t. This flat t dependence is probably related to the fact that the 
lattice-QCD results are evaluated from heavy quarks (corresponding to a pion 
mass of 352 MeV), the extrapolation to the chiral limit is not yet possible in 
general. Similar observations can also be made for higher moments and for 
the moments of the isoscalar GPD combination. In passing we remark that 
an analogous comparison for the Dirac and Pauli form factors [24] (or higher 
moments of H and E) reveals similar results. It is also to be stressed that, 
with regard to the heavy pion used in present lattice calculations, utmost 
caution is advisable in any comparison between lattice and phenomenological 
results. 
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Figure 2: The axial (left) and the pseudoscalar (right) nucleoli form factors 
versus —t. The solid circles represent the lattice results [26] for a pion mass of 
352 MeV. The experimental data are taken from [25] (presented as a dipole 
fits with the errors shown as the green band) and from [27] (open triangles). 
The form factors evaluated from our GPDs are shown as solid lines. The 
dashed and dotted lines represent the pion-pole (with A at = 0.44 GeV) and 
non-pole contributions to Fp, respectively. (Color online) 



3.2 E 

As already mentioned this GPD consists of two parts, a contribution from 
the meson-pole and a non-pole term; the latter has been ignored in most 
studies (cf. for instance [8, 9, 10, 11]). The pole term can be written as [10] 

El% = 6(| x |< $p {^fj (20) 

where Fp°g is the pole contribution to the pseudoscalar form factor for the 
p — > B transition 5 . The non-pole contribution, E D _ P _, to E is parameterized 
as in (12) - (14). The parameters are quoted in Tab. 1. The power = 
is taken as 5. 

As a first check of our parameterization we evaluate the pseudoscalar form 
factor for the p — >■ n transition, traditionally denoted by Fp 

F P (t) = F^ ole (t) + j^dxE<£l{x^t) . (21) 

5 The meson-exchange pole also contributes to the transversity GPDs Ht, Et and Et- 
The pole contribution to the amplitudes for transversely polarized photons are non-zero. 
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With the help of PCAC and the Goldberger-Treiman relation the pion-pole 
contribution to the pseudoscalar form factor can be expressed by 



Ff c (t) = 2mf n ^9*pnF*pn(t) (22) 



m% — t 

where f w (— 131 MeV) is the pion decay constant , g npn (= 13.1) the pion- 
nucleon coupling constant [28] and 

_ A-N — m l / 9 o\ 
r npn — ^2 _ £ ■ \^°> 

The results for Fp(t) are shown in Fig. 2 and compared to the experimental 
data [27] and results from lattice QCD [26]. One notes the same features 
as for the axial form factor - reasonable agreement at low —t and a flatter t 
dependence of the lattice results. The dominance of the pion pole at small t 
is clearly visible. 

From the analysis of ir + electroproduction [1] only the isovector combi- 
nation of i?n.p. is fixed. In the absence of any experimental information that 
would allow to fix the normalization of the isoscalar combination too we have 
to rely on lattice-QCD results. Adjusting the normalizations to the second 
moments of E given in [26] which are free from the pion-pole contribution, 
we arrive at the values of the normalization constants for individual flavors 
quoted in Tab. 1. The lattice results [26] are however not accurate enough to 
allow for a determination of further differences in the parameterization (13) 
of E n . v . for u and d quarks and, we repeat, are not extrapolated to the chiral 
limit. 



3.3 H r £ 

The forward limit of H?, the transversity PDF, is taken from an updated 
version of the analysis of the data on the azimuthal asymmetry in semi- 
inclusive deep inelastic scattering and in inclusive two-hadron production 
in electron-positron annihilation [29]. Since the parameters quoted in [29] 
differ from those obtained in a previous anlaysis which has been used by us 
in [1], our parameterization of H T is to be changed accordingly. Now, we 
parameterize the transversity PDFs as 

5 a = N a HT p 1 / 2 (l-p)[ q a (p) + Aq a (p)]. (24) 
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This parameterization corresponds to a Regge intercept of about zero. The 
PDFs are taken from [23, 30]. The GPD is then calculated from (12) and 
(14) with the Regge parameters quoted in Tab. 1. 

Moments of Ht have been calculated in [31]. They are represented as 
dipole fits with extrapolations to the physical pion mass. In Tab. 2 the 
two lowest moments at t — are compiled and their t dependence is shown 
and compared to the moments evaluated from our GPDs in Fig. 3. Our 
moments at t — which are practically those obtained by Anselmino et al 
[29], are about 40% smaller than the lattice results. As for the other GPDs 
our moments possess a somewhat steeper gradient than the lattice moments 
as, at least for u quarks, can be seen from Fig. 3. The first moment of H T , 
termed the tensor form factor, 

<#(*) = J^dxH^U), (25) 

has also been calculated within a chiral quark-soliton model [32]. At t — 
the following values have been obtained in [32] 

i££ 10 (0) = 0.876, #tio(°) = -0.251, (26) 

which are close the lattice results, and therefore larger than our ones, too. On 
the other hand, the t dependence obtained within that chiral quark-soliton 
model is even steeper than we found, see Fig. 3. Several other models (cf. 
[29] and references therein) also predict larger moments of the transversity 
PDFs than those obtained from the parameterization (24) which is based on 
the analysis performed in [29]. 



3.4 E T 

The GPD Et defined in (1), is a new element in our analysis, it has been 
ignored in [1]. According to lattice QCD [14] Et has however large moments 
with the same sign and similar size for u and d quarks. It therefore seems 
to be expedient to examine its impact on electroproduction of pseudoscalar 
mesons. 

The GPD Et is parameterized as the other ones, (12) - (14), but in 
contrast to the situation for the other GPDs, its parameters can only be 
fixed with the help of lattice results. The latter are presented in [14] in form 
of a p-pole fit 

Tno{ ) ~ [i - t/(pmi n n (27) 
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[31] 


this work 




[14] 


this work 


TTU 

-"no 


0.857(13) 


0.585 


^TIO 


2.93(13) 


2.93 


TTU 

rl T20 


0.268(6) 


0.123 


J2j T20 


0.420(31) 


0.360 


TTd 

n no 


-0.212(5) 


-0.153 


fpd 


1.90(9) 


1.90 


TTd 

rL T20 


-0.052(2) 


-0.021 


fpd 
- Cj T20 


0.260(23) 


0.199 



Table 2: The first two moments of the transversity PDFs Ht and Et at t — 
quoted at the scale 2 GeV. 

with p = 2.5. The forward values, quoted in Tab. 2, are extrapolated to 
the chiral limit but not the mass parameters. The lowest moments at t — 
are the analogues of the familiar anomalous magnetic moments, K, a , and 
are therefore termed tensor anomalous magnetic moments, k t - The lattice 
results for them are of the same magnitude as the usual ones but in contrast 
to those they are both positive. This is in agreement with model studies 
performed in [33] which also support the expectation that k t > | K a \ . The 
lattice result is further corroborated by large-iV c considerations [33] which 
provide k t ~ k t and by a recent constituent quark model [34]. It has been 
speculated [33] that Et is linearly related to the Boer-Mulders function, hj; a , 
[35] while the Sivers function, /" T , is analogously related to the transversity 
PDFs. In a recent analysis [36] of the (cos (20)) measurements [37, 38] some 
evidence has been obtained that the signs of the Boer-Mulders and the Sivers 
functions exhibit the same pattern as those of the lowest moments of Et and 
H T at t — which are quoted in Tab. 2. 

The parameters of E T are fixed by us in such a way that the lowest mo- 
ments from lattice QCD at t — are exactly reproduced and the second 
moments at least approximately. The t dependence is assumed to be some- 
what steeper than the lattice results as is the case for the other GPDs. This 
can be achieved with the parameters 

0k =4, /4 = 5. (28) 

and those quoted in Tab. 1. The moments at t — are given in Tab. 2 and 
their t dependencies are shown in Fig. 3. 
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Figure 3: The first and second moments of (left) and (right) for u 
quarks scaled by the forward value. The results evaluated from our GPDs 
are shown as solid and dashed lines, the lattice-QCD results [14, 31] as bands 
indicating their uncertainties. The dotted line represents the chiral quark- 
soliton model [32]. (Color online) 



We have checked that our GPDs respect various positivity bounds. Thus, 
for instance, 

(E±E T ) 2 (t = 0) <4m 2 (H -H)^(H + H±2H T )(t = 0), (29) 

and similar ones derived in [39] at zero skewness or bounds on H and E in 
terms of PDFs [40] as for instance 



2m 



< - 
~ 4 



(q + Aq + 25) x (q - Aq), 



+ y/{q-Aq) s (q + Aq + 2S) x > 
+ \/{q + Aq-25) x {q- Aq) x , 



+ J(q-Aq) x (q + Aq-26) 3 



(30) 



where the GPDs are to be taken at x, £ and t while the PDFs are either to 
be evaluated at i or i', the individual momentum fractions of the emitted 
or reabsorbed partons, 

x + £ , x — £ 



x 



X 



(31) 
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0.6 
0.4 

ts. 2 

i 0.0 

w i- 

< -0.2 
-0.4 

"°-%.0 0.2 0.4 0.6 

-t'[GeV 2 ] 

Figure 4: Results for the sin (0 — S ) moment for 7r + production shown as 
solid line. The dash-dot-dotted line is obtained under neglect of Ex- Data 
are taken from [6]. (Color online) 



All bounds hold for a given flavor; (30) only holds in the DGLAP region. 
Flavor labels are omitted for convenience in (30). 

4 Electroproduction of pions 

Since we have slightly changed the parameterizations of the GPDs as com- 
pared to [1] and have taken into account now, we have to repeat the 
computation of tt + electroproduction in order to examine the quality of the 
new results. As a matter of fact, it turns out that they are very close to the 
previous results presented in [1]. Their agreement with experiment [5, 6, 7] 
is reasonable. The contribution from E? to the cross section is tiny, much 
smaller than the errors of the HERMES cross section data [5] , since the flavor 
combination Ej, — E^, occuring in 7r + electroproduction, is small. It is only 
noticeable in the asymmetry parameters where it improves the agreement 
with experiment in general. As an example the sin (0 — S ) moment of the 
cross section measured with a transversely polarized target [6] is shown in 
Fig. 4 (0 is the azimuthal angle between the lepton and the hadron plane 
while (f) s specifies the orientation of the target spin vector with respect to the 
lepton plane). 

As a second example the unseparated cross section as well as its longi- 
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Figure 5: Results for the 7r + (left) and 7r° (right) cross sections. The solid 
(dashed, dash-dotted) lines represent the unseparated (longitudinal, trans- 
verse) cross sections. The n + data are taken from [5]. (Color online) 



tudinal and transverse components (da = dax + sdaif) for 7r + production 
are shown in Fig. 5 at Q 2 = 3.44 GeV 2 and W = 3.83 GeV. For the ratio 
of the longitudinal and transverse photon fluxes, e, a value of 0.8 is taken. 
Evidently, the agreement with the HERMES data [5] is very good. The lon- 
gitudinal cross section is dominant at low —if but if —if becomes larger than 
about 0.2 GeV the transverse cross section takes the lead. A crossing of 
the two cross sections is also seen in the large-skewness Jefferson Lab data 
measured by the F n collaboration [13] although it occurs at a smaller value 
of —t'. Actually the crossing takes place at about the same value of t for 
our results and the F n data. It is to be stressed that the dominance of the 
longitudinal cross section at small —t' is a consequence of the strong pion 
pole. Without it one would have dax > da^ at all if. The pure pole contri- 
bution to the longitudinal cross section behaves as —tj (t — m 2 ) 2 at small —t. 
This factor has a pronounced maximum at t — —m 2 . For £ > m n /(2m) the 
maximum lies outside the scattering region. Hence, the longitudinal cross 
section decreases continuously. On the other hand, for £ < m 7r /(2m) the 
above factor generates a maximum of da^/dt at a small value of —if. The 
maximum is visible in the cross section if Q 2 < (m 7r /m)W 2 . 

Before we turn to the discussion of 7r° production it is in order to com- 
ment on the pion wave functions which are needed in the calculation of the 
subprocess amplitudes (cf. (7)). In fact we use the same wave functions as 
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in [1], namely a Gaussian in kj_ with an associated distribution amplitudes 
&as — 6r(l — r) for the twist-2 subprocess amplitudes and $ p = 1 for twist 
3. For the transverse-size parameter in the Gaussian we use a n = [v^tt/V] -1 
(= 0.861 GeV" 1 ) for twist 2 and a wp = 1.8 GeV" 1 for twist 3. There 
is an extra factor multiplying the Gaussian for the twist-3 wave func- 
tion. For a detailed discussion of such a factor and its interpretation see 
[41, 42]. Higher-order Gegenbauer terms in the distribution amplitudes are 
strongly suppressed in the modified perturbative approach: The Sudakov fac- 
tor in conjunction with the subprocess amplitude provides a series of power 
suppressed terms which are generated in the region of soft quark momenta 
(t, 1 — t — y 1). They grow with the Gegenbauer index and reduce the per- 
turbative contribution. With increasing Q 2 the higher Gegenbauer terms 
become gradually more important (cf. the discussion in [43]). Since we are 
merely interested in rather small values of Q 2 the asymptotic distribution 
amplitude suffices for pion production. 

Now, for 7T° electroproduction there is no contribution from the pion pole 
and the GPDs appear in the flavor combination 

(e u F? - e d F?)/V2 (32) 

in contrast to Ff — Ff for tt + production 6 . From the relative signs and sizes 
of the GPDs for u and d quarks it is evident that the contributions from H 
and Ht are large for n + but small for 7r° production while for E mp , and Et the 
situation is reversed. The contributions from the latter two GPDs are small 
for 7r + but large for 7r° production. As a consequence the longitudinal n° 
cross section is much smaller than the corresponding ir + one even if the pion- 
pole contribution is subtracted. Another consequence is that 7r° production 
is dominated by contributions from transversely polarized photons which are 
mainly generated by Et- The latter contribution is by order of magnitude 
stronger in 7r° than in 7r + production where a strong cancellation between E T 
and E T occurs. The importance of E? in 7r° production is obvious from the 
forward dip in the cross section (see Fig. 5); it only contributes to helicity 
flip amplitudes, see (11). 

As we repeatedly mentioned our GPDs are optimized for the range of 
small skewness. Therefore, our approach cannot readily be applied at the 
kinematics accessible at the Jefferson Lab. Ignoring this supposition and ap- 

6 For the production of n + and kaons the quark charges are absorbed in the hard 
scattering amplitude. 
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Figure 6: Unseparated and partial cross sections for tt° production. The 
solid (dashed, dash-dotted) line represents the unseparated (longitudinal- 
transverse, transverse-transverse) cross section. The dash-dot-dotted line is 
evaluated from the modified GPDs, see text. (Color online) 



plying our approach in its present form at Jefferson Lab kinematics straight- 
forwardly, one runs into difficulties with p° production - the cross section is 
at variance with experiment [44] by order of magnitude [2, 45]. For produc- 
tion [46], on the other hand, our approach works quite well. The reason for 
this failure of the handbag approach with p° production is not understood 
but it is likely be related to the valence quarks. Similar difficulties may hap- 
pen for electroproduction of pions. With this admonition in mind one may 
compare our approach with the preliminary CLAS data on ir° production 
[47] (at W ~ 2.5 GeV and Q 2 ~ 2.3 GeV 2 ). It turns out that the magnitude 
of the cross section is about right. Even a little forward dip is seen in the 
CLAS data on the cross section. This effect has also been observed in a 
Hall A experiment [48]. It is important to realize that the depth of the dip 
is determined by the ratio of the amplitudes .Mo-,++ an d Ai 0+t++ , i.e. it is 
influenced by the relative strength of and Et- We remind the reader that 
the normalization of the first GPD is fixed by the transversity PDF deter- 
mined in [29] while the second one is constrained by lattice-QCD results [14]. 
A change of these GPDs may lead to a different shape of the cross section in 
the forward region. An example of such a modification of the GPDs is also 
shown in Fig. 6 and compared to the result evaluated from the standard pa- 
rameterization of the GPDs presented in Sect. 3. For this variant we modify 
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Figure 7: Results for 7r° the cross section versus t' at W = 5 GeV for sample 
values of Q 2 (left) and at Q 2 = 3 GeV 2 for various values of W (right). The 
polarization of the photon flux is taken as e = 0.35 at W = 3.5 GeV and 0.8 
for W > 5 GeV. (Color online) 



Ht by changing its normalization: Njj = 1.1 and Nfj T = —0.3 and take 
for the slope bn T = 0.3 GeV~ 2 (see Tab. 1). In order to have more or less 
the same 7r° cross section at CLAS kinematics the GPD Ej, is multiplied by 
a factor 0.8. For this modification of the transversity GPDs our fit to 7r + 
production cross section practically remains unchanged. It is of the same 
quality as for the standard parameterization since — is almost unal- 
tered. On the other hand, the forward limits of Ht are now at variance with 
the transversity PDFs given in [29] and the lowest moments of are below 
the lattice-QCD results [14]. As one sees from Fig. 6 the dip at t' is consid- 
erably less deep now. The longitudinal-transverse and transverse-transverse 
interference cross sections are also shown for the standard parameterization 
in this figure. The analogous predictions for CLAS kinematics are in fair 
agreement with experiment [47]. The energy and Q 2 dependencies of the ir° 
cross section are shown in Figs. 7. 

Predictions for the 7r° target asymmetries are shown in Fig. 8 for the stan- 
dard kinematics as an example. The transverse target asymmetries behave 
differently for 7r° and 7r + production. Particularly noteworthy is the fact the 
sin0 s moment for 7r° production is much smaller than for tt + for non-zero 
t'. This behavior can be traced back to the amplitude A^o+,++ which is very 
small for tt + but large for 7r° production. In terms of helicity amplitudes the 
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Figure 8: Results for the sin (0 — <p s ) and smcj) s moments of the transverse 
target asymmetries for 7r° electroproduction versus t' . (Color online) 



asymmetry reads 



A^a ~j£(l+e)ImMt 



M 



o-,o+ 



Ml ++ M 



o+,o+ 



(33) 



where it is assumed that the angle of the rotation in the lepton plane from 
the direction of the incoming lepton to the virtual photon is small. The 
quantity do is given by 



<TQ 



1 r 
2 



+ 



-^0+,+ 

M 0+ ^ 



M -,- 



Mo-, 



+ I Mq+.-a 



+ I M 



o-,o+ 



(34) 



The 7r + asymmetry is essentially generated by the second term in (33) while, 
for not too small —t', a strong cancellation between both the terms in (33) 
occur for ir° production. For t' — > the first term vanishes oc t'. 

The different behavior of the sin (<p — <p s ) moment for ir + and tt° pro- 
duction (see Figs. 4 and 8) is also a consequence of the contribution from 
Et being small in one but large in the other reaction. Approximately this 
asymmetry is given by 



* sin (</>-(/>„) 



-2elm 



Im 



Ml, 



.Mo-, 



(35) 



The second term provides a positive contribution to the ir° asymmetry which 
overcompensates the asymptotically dominant first term provided by longi- 
tudinally polarized photons. 
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Our estimates of ir° electroproduction which, we repeat, strongly depend 
on the lattice QCD findings [14] for the moments of Et, are distinct from 
our previous estimate presented in [1] where the amplitude A^o+,++ had been 
neglected, and from other predictions. In collinear factorization to leading- 
twist accuracy only the longitudinal cross section can be calculated. The 
predictions for it, see e.g. [10], are of about the same size as our longitu- 
dinal cross section but much smaller than the full cross section. In [49] the 
dominance of the transverse cross section in it production has also been sug- 
gested. However, our treatment of the twist-3 contribution differs from the 
approach advocated in [49] markedly. In the latter work the subprocess is 
viewed as form factors for photon-pion transitions under the action of vector 
and axial-vector currents. A proportionality between E? and Ht is assumed 
where the constants of proportionality are set by the transverse anomalous 
magnetic moments for which the values 0.6 for u and 0.3 for d quarks are 
utilized. This parameterization of E T is not supported by lattice QCD [14] 
(see Tab. 2) and differs from our parameterization drastically: E T in [49] is 
much smaller than our one and has opposite signs for u and d quarks. As 
a consequence the absolute value of the amplitude Ai _ i++ is much larger 
than that of A^o+,++ as opposed to our results. The results for the 7r° cross 
section are therefore substantially smaller in [49] . Instead of a forward dip a 
forward maximum occurs. 



5 Electroproduction of of 77 and rj 

The treatment of the 77 and rf mesons within the handbag approach is similar 
to the case of the 7r°. It is only a bit more intricate due to rj — rj' mixing and, 
in principle, the gluon-gluon Fock component of the rj and rj' mesons. This 
Fock component contributes to the same order of a s as the quark- ant iquark 
components as can be seen from the Feynman graph shown in Fig. 9. It 
has however been shown in [50] that the twist- 2 subprocess amplitude for 
the gluon-gluon component is oc t/Q 2 and is therefore to be neglected for 
consistency. A possible twist-3 gluon contribution is neglected as well since, 
in contrast to the quark subprocess (cf. (9)), there is no known reason why 
this contribution could be large. 

Working in the flavor octet-singlet basis and exploiting the results of the 
quark-flavor mixing scheme developed in [51], we can decompose any of the 
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Figure 9: A typical Feynman graph for the contribution from the gluon-gluon 
Fock component of the r\ and rj' mesons. 



helicity amplitudes as 

M v = cos6 8 M (8) -sin0iM (1) , 
M v> = sin# 8 .M (8) + cosfl 1 .M (1) . (36) 

For the mixing angles we adopt the values 

# 8 = -21.2°, 6 1 = -9.2° (37) 

derived in [51] on exploiting the divergencies of the axial vector current which 
embodies the axial-vector anomaly. The octet and singlet amplitudes are to 
be calculated in full analogy to the case of the tc° with octet and singlet wave 
functions and decay constants taken from [51] 

/ 8 = 1.26/ w , /i = 1.17/ ff . (38) 

For these wave functions, assumed to be independent on the meson, we take 
the same Gaussian as for the pion with exactly the same transverse size pa- 
rameters. For Q 2 in the range of interest in this work these simple wave 
functions suffice, cf. the remark in Sect. 4. In the analysis of the Babar data 
on the P'-f transition form factors [52], also performed within the modified 
perturbative approach [43], only little differences between the octet and sin- 
glet wave functions have been found. The GPDs appear in the following 
flavor combinations 

FP = 4= [zuF? + e d F? - 2e s F-} , 

FP = [eJF? + e d Ff + e.Ff\ . (39) 
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Figure 10: Results for the 77 (left) and rj' (right) cross sections. For notations 
refer to Fig. 5. (Color online) 



For the charge-conjugation-even mesons the GPDs contribute in the valence- 
quark combination F" — F". For strange quarks we assume F* ~ F/ for all 
GPDs 7 . Hence, there is no contribution from strange quarks and one arrives 
at the relation f[ 1] = y2F/ 8) . 

The 7] and rf cross sections are shown in Fig. 10. In order to facilitate 
comparison of the cross sections for the various mesons we use the same 
kinematical setting, W = 3.83 GeV and Q 2 = 3.44 GeV , for which there 
are tt + data from HERMES [5]. For an understanding of the differences 
between the 77 and 77° cross sections it is important to realize that the u and 
<i-quark GPDs contribute with opposite signs in these two reactions (see (39) 
and (32)). For small —t', the region where the GPDs H and H? dominate, 
the opposite sign of their u and <i-quark parts leads to a larger contribution 
to rj than to 77° production which, in collaboration with the larger decay 
constant (38), overcompensates the relative factor l/y/3 between (39) and 
(32). Hence, the 77/77° ratio of the cross sections is of order 1, see Fig. 11. For 
large —t' the GPD Et dominates which has the same sign for u and d quarks. 
Therefore, the 77/77° ratio is much smaller than 1; in fact close to 1/3 for a 
large range of if . Interestingly the CLAS collaboration [47] has measured 
the 77/77° cross section ratio for —t reaching from 0.14 to 0.77 GeV 2 but for 

7 For the unpolarized strange quark PDFs [53] and for H s [54] possible differences 
between strange and antistrange distributions have been studied. No evidence for a non- 
zero difference has been found within errors. 
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Figure 11: Ratio of the rj and n° cross sections. Solid (dash-dot-dotted) line 
represents the results evaluated from the standard parameterization of the 
GPDs (modified Ht). For notations refer to Fig. 5. (Color online) 



values of Q 2 and W that are somewhat too low for allowing a comparison 
with our results without reservation. Nevertheless the preliminary CLAS 
data match very well our predictions which can be regarded as an indication 
of large contributions from the GPD Et with the same sign for its u and 
d-quark parts. We note in passing that our i]/it ratio of the longitudinal 
cross sections at low —t' is in agreement with an estimate presented in [55]. 

The transverse cross section for rj' production is much smaller than those 
for the r) and 7r° channels. The reason is obvious. The twist-3 mechanism is 
not enhanced by the chiral condensate for the flavor-singlet part. Actually 
it is about half as strong as for the octet channels. 

In Fig. 12 the sin (0 — <p s ) and the sin0 s moments of the 77 and 7/ cross 
sections measured with a transversely polarized target are shown. The trends 
of these asymmetries for 77 production bear similarities to the corresponding 
7T° ones while for the 7/ the suppression of the twist-3 effect reflects itself in 
a different behavior of the asymmetries. 

6 The kaon-hyperon channels 

The analysis of kaon electroproduction is similar to n + production; the same 
expressions hold for the convolutions (e<j = e s ). To describe the if-meson we 
again use a Gaussian wave function with a transverse-size parameter ax = cl-k 
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Figure 12: Results for the sin (<ft — <p s ) and sm<f) s moments of the transverse 
target asymmetries for rj (left) and 7/ (right) electroproduction versus t'. 
(Color online) 



and associated distribution amplitudes, the flat one at twist-3 level [56] and 

$ k (t, fi) = $as(t) [l + a K1 (fi)C 3 1 /2 (2r - 1) + a K2 (fi)C 3 2 /2 (2r - 1) + . . . ] (40) 

for twist 2 (the momentum fraction r refers to the u quark). The impor- 
tant difference to the case of pions is that this distribution amplitude is not 
symmetric in the momentum fractions r and 1 — r in general. A restriction 
to the asymptotic distribution amplitude is inappropriate for kaons since it 
would lead to the relation 

U K = f K /UU* (41) 

among the twist-2 subprocess amplitudes. With f K = 1.21f n the kaon am- 
plitudes would then be larger than the pion ones. This is against experience; 
kaon channels are typically suppressed by about 10% against pion channels. 
This can be seen, for instance, in the electromagnetic form factors [57], in 
two-photon annihilations [58] or Xcj decays into pairs of kaons or pions [59]. 
For the lowest two Gegenbauer coefficients we use the values 



a K1 (l GeV) = -0.1, a K2 (l GeV) = -0.1, (42) 

which are in agreement with the Xcj — > KK decay width (J = 0,2) [60]. 
QCD sum rules [56], on the other hand, provide a K i = —0.06 ± 0.03 and 
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0'K2 = 0.30 ± 0.15. The negative values of clki imply that, as expected, 
the antistrange quark in the kaon carries a larger momentum fraction on 
the average than the non-strange one ((r) = 1/2 + 3/lOa^i; (1 — r) = 
1/2 — 3/IOor-i). We take the first two Gegenbauer terms with the coefficients 
(42), into account in order to avoid an overestimate of kaon production. Their 
suppression within the modified perturbative approach (cf. the discussion in 
Sect. 4) is mild. For the twist-3 distribution amplitude we ignore such an 
effect. 

For kaon production the proton-hyperon transition GPDs enter which, 
with the help of SU(3) flavor symmetry, can be related to the proton GPDs 
by [9] 

Fip^-A — 

F ip ^+ ~ 

As for the other channels we assume a flavor-symmetric sea for all GPDs 
and Ff — F? ~ 0. In contrast to the case of the 7r + where the transition 
GPDs Fi P ^ n are related to the proton ones by isospin invariance, one has to 
be aware of possible large flavor-symmetry breaking effects in (43). 

Finally, we have to discuss the kaon pole. It is treated in analogy to the 
pion pole (see left hand graph in Fig. 1). The contribution of the kaon pole 
to Ep^y is given in (20) and that to the p — > Y form factor reads 

where the form factor Fk p y parameterizes a residual t dependence like the 
one which occurs for the pion. For this form factor which has not been 
considered in previous work (see for instance [61]), one may take the same 
expression as for the case of the pion (23) for a first estimate. The GPD E 
satisfies the sum rule [61] 

\[ x dxE p ^y(x, £, t) = F KY (t) + -^ Y (t) . (45) 

Due to violation of flavor symmetry by the strange quark mass the extra form 
factor Y occurs in the sum rule. This of no relevance here since the sum 



± i ~ ± i i 



Ff 



(43) 
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rule (45) will not be exploited. It is important to realize that in (20) and (45) 
the transition GPD actually appears, no assumption on flavor symmetry is 
required. For the kaon-baryon coupling constants appearing in (44) we use 
values obtained from fixed-t dispersion relations which are however in good 
agreement with SU(3) predictions [62] 

9k+ p a ^ -13.3 , g K + P T,° = -9k» p y.+ /^ ^ -3.5 . (46) 

The relation between the coupling constants for K + pYP and K G pL + is a con- 
sequence of isospin invariance. The residue of the kaon pole (5) is fully fixed 
by noting that we take for the electromagnetic form factor of the positively 
charged kaon that of the pion but multiplied by a factor of 0.9 in order to 
take care of flavor symmetry breaking 

9 

Fk+ m = 1 + Q2/Q.462 GeV 2 " (4?) 

This is at least in agreement with the CLEO measurement [57] of the time- 
like pion and kaon electromagnetic form factors. 

The electromagnetic form factor of the K° is not necessarily zero. In 
the perturbative frame work for instance it is related to the antisymmetric 
Gegenbauer coefficients in (40). Since nothing is known about this form fac- 
tor we assume its ratio to the K + form factor to be given by the perturbative 
ratio 

3 1 + a 2 K1 + a K2 + 2a K2 

Evaluating this ratio from the Gegenbauer coefficient quoted in (42), one 
finds a value of -0.15. 

Our results for cross sections and target asymmetries are shown in Figs. 
13 and 14. The kaon cross sections are dominated by contributions from 
transversely polarized photons too. As opposed to n + production the meson- 
pole term plays a minor role here because of the large kaon mass appearing 
in the pole denominator. It is however strong enough to produce a little 
maximum in the longitudinal cross section for K + A at t' ~ —0.1 GeV 2 by 
means of the factor —t/(t — m 2 K ) 2 in the pure pole contribution {m 2 K > —to)- 
For the two other kaon channels the pole contribution is hardly visible. The 
cross section for the K + A channel is substantially larger than those for the 
other kaon channels. Besides different pole contributions, these two processes 
differ from each other by the factor 1 / v^2 in the GPD combinations (43) in 
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Figure 13: Results for the cross sectios of the K + A (left), K + TP (center) and 
iT°E + (right) final states. For notations refer to Fig. 5. (Color online) 




Figure 14: Results for the sin (<ft — <p s ) and sin <p s moments of the transverse 
target asymmetries for electroproduction of the K + A (left), K + TP (center) 
and (right) final states versus if . (Color online) 



favor of 7*p — > K°Ti + . This factor is overcompensated by the w-quark charge 
in the hard amplitude of the process 7*p — > K + T! ) instead of in the other 
process. As a result the K + TP cross section is about twice as large as the 
one. On the other hand, the transverse target asymmetries behave 
very similar for both the processes, see Fig. 14. The asymmetries for the 
K + A channel are similar to those for 7r + production. 

Finally, in Fig. 15, we compare the various cross sections at a larger value 
of W, characteristic of the COMPASS experiment. 
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Figure 15: The cross sections for various pseudoscalar meson channels at 
W = 7 GeV and Q 2 = 3 GeV 2 . (Color online) 



7 The longitudinal beam and target 
asymmetries 

In Fig. 16 we display the beam-spin asymmetry for the pseudoscalar-meson 
channels at the standard kinematical setting. In terms of helicity amplitudes 
this observable is expressed by 



A LU a Q = ^e(l-e)lm[(M* 0+t++ -M* 0+ ^ + )M 0+ ,o + 

+ (M* _ >++ -M* _ j _ + )M „,o + }, (49) 

In our handbag approach the first term in (49) cancels (see (11)) and the 
second one simplifies to 



A LU a ~ y/e(l - e)lm[M* _ ++ -M -,o+J • (50) 

For reactions being dominated by as for example it production, A^u 
is very small. Others like 7r + production for which E? is less prominent, 
are large. This can be clearly seen in Fig. 16. The CLAS collaboration 
[63] has measured Am for tc° and n + production at W ^ 2.4 GeV and 
Q 2 ~ 2.5 GeV 2 . While the CLAS result for % + production is in reasonable 
agreement with our result, is that for tt° production larger than our one. 

Another interesting observable is the sin 0-moment of the electroproduc- 
tion cross section measured with a longitudinally polarized target. To the 
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Figure 16: The beam spin asymmetry (left) and the asymmetry for a lon- 
gitudinally polarized target (right) for various pseudoscalar-meson channels 
versus if. (Color online) 



extend that the angle of the rotation in the lepton plane from the direction 
of the incoming lepton to the virtual photon is small the longitudinal target 
asymmetry reads 



A UL a ~ -^e(l + e)lm[(M* 0+)++ + M* 0+r 
+ (A4S-.++ + M *o-,-+) M)-,o+ 
It simplifies in our handbag approach to 

A ul <tq ~ 



M 



0+,0+ 



(51) 



e(l + e)Im 



2M1 



M 



0+,(H 



+ M 



0-,++ 



M - 0+ 



(52) 



Without the GPD the ratio of the two asymmetries would be given by 
Aul/Alu = 



(1 + s)/(l — e). Deviations from this ratio signal a non-zero 
contributions from Et, see Fig. 16. 



8 Summary 

We have investigated electroproduction of pseudoscalar mesons within the 
handbag approach restricting ourselves to the range of small skewness. Forced 
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by the behavior of target asymmetries for 7r + production measured by the 
HERMES collaboration [6] we have allowed for a twist-3 effect consisting of a 
twist-3 meson distribution amplitude in combination with the leading-twist 
transversity GPDs. This twist-3 effect is strongly enhanced by the chiral 
condensate (9), (10). It therefore seems to be justified to include this effect 
in the analysis even though we don't perform a systematic analysis at twist-3 
accuracy. We stress that for electroproduction of vector mesons the twist-3 
effect is expected to be much smaller since /ip is to be replaced by the vector 
meson's mass. 

The GPDs are constructed from double distributions which are parame- 
terized as zero-skewness GPDs multiplied by a Regge-like t dependence and 
a weight function from which the skewness dependence of the GPD is gen- 
erated. They are constrained by form factors and PDFs. Their isovector 
combination are fitted to the HERMES data on 7r + electroproduction. Mo- 
ments of the GPDs are compared to moments obtained from lattice QCD. 
In general reasonable agreement is to be observed. A number of positivity 
bounds for the GPDs are checked as well by us. An exceptional case is Et- Its 
parameterization is fully constrained by lattice-QCD results [14] since this is 
the only available information about it at present. Therefore, Et suffers from 
the usual uncertainties of lattice-QCD results arising from the still low statis- 
tics and from the compulsion of working with heavy quarks. It turns out that 
the GPD Et for u and d quarks cancel to a high degree in n + electroproduc- 
tion while it plays a prominent role in most of the other pseudoscalar-meson 
channels. In fact, with the exception of tt + and rj' they are all dominated 
by contributions from the transversity GPDs, i.e. by transversely polarized 
photons, in the range of Q 2 accessible in present-day experiments. Of course, 
for Q 2 — > oo the longitudinal cross sections will dominate in all reactions. 
These findings are to be contrasted with the usual assumption made in the 
handbag approach [9]- [12] that contributions from longitudinally polarized 
photons already dominate for Q 2 of the order of a few GeV. 

A precise calculation of all details of electroproduction of pseudoscalar 
mesons is beyond feasibility at present. There are many uncertainties like 
the parameterization of the transversity GPDs or the exact treatment of the 
twist-3 contribution (e.g. the neglect of possible three-particle configuration 
of the meson state). Also higher-order perturbative corrections other than 
those included in the Sudakov factor and, implicitly, in the experimental 
electromagnetic form factors of the mesons, are ignored. According to [64] 
the NLO corrections are rather large for the cross sections for Q 2 < 10 GeV. 
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However, the results presented in [64] refer to the leading-twist contribution 
to 7r + production, i.e. to about 10% of the cross section since only the per- 
turbative contribution to the pion form factor is taken into account. Thus, 
it is not clear how strong the NLO corrections are in our approach. Further 
uncertainties occur for kaon production. In contrast to the case of the n + 
where the p — > n transition GPDs are related to the diagonal proton ones 
by isospin symmetry, the proton - hyperon transition GPDs are connected 
to the proton GPDs by SU(3) flavor symmetry which is less accurate than 
isospin symmetry. The assumption of a flavor symmetric sea for all GPDs is 
also stronger for kaons than for pions. With regard to all these uncertainties 
we consider our investigation of electroproduction of pseudoscalar mesons 
as an estimate of the pertinent observables. The trends and magnitudes of 
the cross sections are likely correct but probably not the details. Particularly 
large uncertainties are to be expected for subtle observables as the spin asym- 
metries. On the other hand, the large magnitude of contributions from the 
transversity GPDs seems to us hard to avoid. Future data measured at the 
upgraded Jefferson Lab or by the COMPASS collaboration may probe our 
results. Although our approach is designed for small skewness (and rather 
large Q 2 and W) we also compared our results with the large skewness data 
measured at the present JLab accelerator. We found general agreement in 
trends and magnitudes. We regard this as a hint at strong contributions from 
the transversity GPD E T for which the u and d quark parts have the same 
sign. 
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